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1 Introduction 



In the imminent LHC environment, where one expects to have an experi- 
mental luminosity precision tag at the level of 2%, [1] the requirement for 
the theoretical precision tag on the corresponding luminosity processes, such 
as single W, Z production with the subsequent decay into light lepton pairs, 
should be at the 0.67% level in order not to compromise, unnecessarily, the 
over-all precision of the respective LHC luminosity determinations. This dic- 
tates that multiple gluon and photon radiative effects must be controlled at 
the stated precision. The theory of QED ® QCD exponentiation [2] allows 
for the simultaneous resummation of multiple gluon and multiple photon ra- 
diative effects in LHC physics processes, to be realized ultimately by MC 
methods on an event-by-event basis in the presence of parton showers, in 
a framework which allows us to systematically improve the accuracy of the 
calculations without double-counting of effects, in principle to all orders in 
both a s and a. Such a theoretical framework opens the way to the desired 
theoretical precision tag on the LHC luminosity processes. 

Our starting point for the new QED(g>QCD resummation theory [2] is 
the QCD resummation theory presented in Ref. [3]. This resummation is an 
exact rearrangement of the QCD perturbative series based on the N — 1 term 
in the exponent in the formal proof of exponentiation in non-Abelian gauge 
theories in the eikonal approximation, as given in Ref. [4]. This exponential 
is augmented with a sum of residuals which take into account the remaining 
contributions to the perturbative series exactly to all orders in a s We 
therefore have an exact result whereas the resummation theory in Ref. [4] and 
those in Refs. [5-7] are approximate. Recently, an alternative resummation 
theory, the soft-collinear effective theory(SCET) [8], has been developed to 
treat double resummation of soft and collinear effects. Since we have an exact 
re-arrangement of the perturbative series, we could introduce the results from 
Refs. [5-8] into our representation as well. Such introductions will appear 
elsewhere. 

The need for the extension of the QCD resummation theory to QED® 
QCD resummation was already suggested by the results in Refs. [9-14], where 

^df desired, our overall expoential factor can be made to include all of the terms in the 
exponent in Ref. [4], in principle. 
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it was shown that in the evolution of the structure functions the inclusion of 
the QED contributions leads to effects at the level of ~ 0.3%, already almost 
half of the error budget discussed above. We will find similar size effects from 
the threshold region of heavy gauge boson production. All of these must be 
taken into account if one wants ~ 1.0% for the theoretical precision tag. 

The discussion is organized as follows. In Section 2, we review the exten- 
sion of the YFS theory to an exact resummation theory for QCD. Section 
3 presents the further extension to QED®QCD. Section 4 contains the ap- 
plication to heavy gauge boson production with the attendant discussion of 
shower /ME matching. Section 5 contains some concluding remarks. 



2 Extension of YFS Theory to QCD 



We consider a parton-level single heavy boson production process such as 
q + q' -> V + n{g) + X -> U' + n{g) + X, where V = W ± , Z, and £ = 
e, /i, £' = v e , V/j,(e, /i) respectively for V = W + (Z), and £ = v e ,v^ £' — e, /i 
respectively for V = W~. It has been established [3] that the cross section 
may be expressed as 



/ 



n=0 J 3=1 ■> 
a y _ e iy(pi+P2-qi-q2-T,kj)+D QCD (1) 



(27T) 
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where gluon residuals /3 n (ki, . . . ,k n ), defined by Ref. [3], are free of all in- 
frared divergences to all orders in a s (Q). The functions SUM/^QCD) and 
-Dqcd, together with the basic infrared functions -Bqcd> ^qcd> an d ^qcd are 
specified in Ref. [3]. We call attention to the essential compensation between 
the left over genuine non-Abelian IR virtual and real singularities between 
the phase space integrals J dPh f3 n and J dPh /3 n+ i that really allows us to 

isolate j3j and distinguishes QCD from QED, where no such compensation 
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occurs. The result in (CQ) has been realized by Monte Carlo methods [3]. See 
also Refs. [15-17] for exact 0(a 2 s ) and Refs. [18-20] for exact 0(a) results 
on the heavy gauge boson production processes which we discuss here. 

Apparently, we can not emphasize too much the exactness of ([T]). Some 
confusion seems to exist because it does not show explicitly an ordered ex- 
ponential operator for an appropriate ordering prescription, path-ordered, 
time-ordered, etc. The essential point is that, in ([1]), we have evaluated 
the matrix elements of these operators and written the result in terms of 
the over- all exponent shown therein and the residuals j3j. This allows us to 
maintain exactness to all orders in a s . 



3 QED0QCD Resummation Theory 



The new QED <g> QCD theory is obtained by simultaneously exponentiating 
the large IR terms in QCD and the exact IR divergent terms in QED, so that 
we arrive at the new result 

da cxp = e SUM «(QCED) £ / 

n,m=0 ^ ji=l 31 j 2 =l J2 

/ dAy r iv(p±+q±-P2-ci2-Y,k^ -Efc'j 2 )+-PQCED (2) 
J W 

v n (b h -y y \ d3p2 d3q2 

where the new YFS [21,22] residuals, /3 nim (fci, . . . , k n ;k[, . . . , k' m ), with n hard 
gluons and m hard photons, defined in Ref. [2], represent the successive 
application of the YFS expansion first for QCD and subsequently for QED. 

The functions SUMir(QCED), -Dqced are determined from their QCD 
analogs SUM IR (QCD), -Dqcd via the substitutions 

TDnls TDnls I TDnls — TDnls 

-°QCD ~^ -°QCD ~r -°QED — -°QCED5 

TDnls finis I finis — fjnls tn\ 

-°QCD — *■ -°QCD "T -°QED — -°QCED5 \°) 

onls onls I cnls — onls 

°QCD ~ * °QCD "T" °QED — °QCED 
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everywhere in expressions for the latter functions given in Ref . [3] . We stress 
that if desired the exponent corresponding the N th Gatherall exponent for 
N > 1 can be systematically included in the QCD exponents SUMir(QCD), 
-Dqcd if desired, with a corresponding change in the respective residuals 
Pn,m(ki, . . . , k n \ k'i, . . . , k' m ) . The residuals P n ,m(ki, . . . , k n \ k±, . . . , k m ) are 
free of all infrared singularities, and the result in is a representation that 
is exact and that can therefore be used to make contact with parton shower 
MC's without double counting or the unnecessary averaging of effects such as 
the gluon azimuthal angular distribution relative to its parent's momentum 
direction. 

In the respective infrared algebra (QCED) in the average Bjorken x 
values 

x avg {QED) = 7 (QED)/(1 + 7 (QED)), 
WQCD) = 7 (QCD)/(1 + 7 (QCD)), 

where 7(A) = ^^(L s - 1), A = QED, QCD, with C A = Qj,C F , respec- 
tively, for A = QED, QCD and the big log L s , imply that QCD dominant 
corrections happen an order of magnitude earlier than those for QED. This 
means that the leading /? 0)0 -level gives already a good estimate of the size of 
the interplay between the higher order QED and QCD effects which we will 
use to illustrate (jSj) here. 



4 QED(g) QCD Threshold Corrections and 
Shower/ME Matching at the LHC 

The cross section for the processes pp — > V + n( r y) +m(g) +X — > ££' + n'{ r y) + 
m(g) +X, where V,£,£' are the vector-boson / lepton combinations defined 
in Section 3, may be constructed from the parton-level cross section via the 
usual formula (we use the standard notation here [2]) 
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In this section, we will use the result in here with semi-analytical methods 
and structure functions from Ref. [23] to examine the size of QED(g>QCD 
threshold corrections. A Monte Carlo realization will appear elsewhere [24]. 

First, we wish to make contact with the existing literature and stan- 
dard practice for QCD parton showers as realized by HERWIG [25] and/or 
PYTHIA [26]. Eventually, we will also make contact with the new parton 
distribution function evolution MC algorithm in Ref. [27]. We intend to 
combine our exact YFS-style resummation calculus with HERWIG and/or 
PYTHIA by using the latter to generate a parton shower starting from the 
initial (xi,X2) point at factorization scale /i, after this point is provided by 
the {Fi}. This combination of theoretical constructs can be systematically 
improved with exact fully exclusive results order- by-order in a s , where cur- 
rently the state of the art in such a calculation is the work in Ref. [28] which 
accomplishes the combination of an exact 0(a s ) correction with HERWIG, 
where the gluon azimuthal angle is averaged in the combination. 

The issue of this being an exact rearrangement of the QCD and QED 
perturbative series requires some comment. Unlike the threshold resumma- 
tion techniques in Refs. [5-7] , we have a resummation which is valid over the 
entire phase space. Thus, it is readily applicable to an exact treatment of 
the respective phase space in its implementation via MC methods. 

We may illustrate how the combination with PYTHIA/HERWIG may 
proceed as follows. We note that, for example, if we use a quark mass 
m q as our collinear limit regulator, DGLAP [29] evolution of the structure 
functions allows us to factorize all the terms that involve powers of the big log 
L c = In fj, 2 /m q — 1 in such a way that the evolved structure function contains 
the effects of summing the leading big logs L = In /i 2 / /i 2 , where the evolution 
involves initial data at the scale /x - This gives us a result independent of m q 
for m q i 0. In the DGLAP theory, the factorization scale \x represents the 
largest pt of the gluon emission included in the structure function. 

In practice, when we use these structure functions with an exact result 
for the residuals in (121), it means that we must in the residuals omit the 
contributions from gluon radiation at scales below /i. This can be shown to 
amount in most cases to replacing L s = In s/m q — 1 — > L nls = lns//i 2 but 
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in any case it is immediate how to limit the pt in the gluon emissions so 
that we do not double count effects. In other words, we apply the standard 
QCD factorization of mass singularities to the cross section in (T5]) in the 
standard way. We may do it with either the mass regulator for the collinear 
singularities or with dimensional regularization of such singularities. The 
final result should be independent of this regulator and this is something 
that we may use as a cross-check on the results. 

This would in practice mean the following: We first make an event with 
the formula in (j3j) which would produce an initial beam state at (xi,X2) for 
the two hard interacting partons at the factorization scale \x from the struc- 
ture functions {Fj} and a corresponding final state X from the exponenti- 
ated cross section in da exp (xiXjs), where we stress that the latter has had all 
collinear singularities factorized so that it is much more convergent then its 
analog in LEP physics for the electroweak theory for example. The standard 
Les Houches procedure [30] of showering this event (xi, x%, X) would then be 
used, employing backward evolution of the initial partons. If we restrict the 
p T as we have indicated above, there would be no double counting of effects. 
Let us call this px matching of the shower from the backward evolution and 
the matrix elements in the QCED exponentiated cross section. 

It is possible, however, to be more accurate in the use of the exact result 
in ([2]). Just as the residuals (3 n>m {kx, . . . , k n ] k[, . . . , k' m ) are computed order 
by order in perturbation theory from the corresponding exact perturbative 
results by expanding the exponents in (jSJ) and comparing the appropriate 
corresponding coefficients of the respective powers of a n a™, so too can the 
shower formula which is used to generate the backward evolution be expanded 
so that the product of the shower formula's perturbative expansion, the per- 
turbative expansion of the exponents in (TJJ, and the perturbative expansions 
of the residuals can be written as an over-all expansion in powers of a n a™ 
and required to match the respective calculated exact result for given order. 

In this way, new shower subtracted residuals, {^ n ,m{ki, ■ ■ ■ , k n ; k[, . . . , k' m )}, 
are calculated that can be used for the entire gluon px phase space with an 
accuracy of the cross section that should in principle be improved compared 
with the first procedure for shower matching presented above. Both ap- 
proaches are under investigation, where we note that the shower subtracted 

2 Here, we refer to both on-shell and off-shell emitted gluons. 
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residuals have been realized for the exact 0{a) luminosity Bhabha process 
at DAPHNE energies by the authors in Ref. [31]. 



Returning to the general discussion, we compute, with and without QED, 
the ratio r exp = o" exp /o"Born, where we do not use the narrow resonance ap- 
proximation, for we wish to set a paradigm for precision heavy vector boson 
studies. The formula which we use for 0"B O rn is obtained from that in (TJJ by 
substituting da-Q OTTL for da cxp therein, where d&Bom is the respective parton- 
level Born cross section. Specifically, we have from (jTJ the /3 ,o-level result 

/""max 

OexpOia^s) = / dv 7qced ^^"^yfsCtqced) 

Jo (5) 

e 5YFS a Bor n((l - v)x 1 x 2 s) 

where we intend the well-known results for the respective parton-level Born 
cross sections and the value of t> max implied by the experimental cuts under 
study. 

What is new here is the value for the QED(g>QCD exponent 

7qced = (2QJ^ + 2CV^) L nls (6) 

where L n j s = lnxiX2s//x 2 when fi is the factorization scale. The functions 
-Fyfs(7qced) and <5yfs(7qced) are well-known [22] as well: 

g— 7QCED7-E 

Fyf s(iqced) 



i\1 + 7qced) ^ 
W7qced) = ^7qced + (Q^ + CV^) (2C(2) - , 

where ((2) is Riemann's zeta function of argument 2, i.e., vr 2 /6, and 7# is 
Euler's constant, i.e., 0.5772. . . . 

Using these formulas in (jl]) allows us to get the results 



1.1901 , QCED = QCD+QED, LHC 

. 1.1872 , QCD, LHC , . 

r ex P -S 11911 ; QCED = QCD+QED, Tevatron ,N 

1.1879 , QCD, Tevatron. 
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We see that QED is at the level of .3% at both LHC and FNAL. This is stable 
under scale variations [2]. We agree with the results in Refs. [15,16,18-20] on 
both of the respective sizes of the QED and QCD effects. Furthermore, the 
QED effect is similar in size to structure function results found in Refs. [9-13]. 

5 Conclusions 

We have shown that YFS theory (EEX and CEEX), when extended to non- 
Abelian gauge theory, allows simultaneous exponentiation of QED and QCD, 
QED<g>QCD exponentiation. For QED<g)QCD we find that full MC event 
generator realization is possible in a way that combines our calculus with 
HERWIG and PYTHIA in principle. Semi-analytical results for QED (and 
QCD) threshold effects agree with literature on Z production. As QED is 
at the .3% level, it is needed for LHC theory predictions at < 1%. The cor- 
responding analysis of the W production is in progress. We have illustrated 
a firm theoretical basis for the realization of the complete 0(a 2 , aa s , a 2 ) re- 
sults needed for the FNAL/LHC/RHIC/ILC physics and all of the latter are 
in progress. 
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